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ABSTRACT 



We present deep Hubble Space Telescope NICMOS ncar-infrarcd and WFPC2 optical imaging of a 
small region in the core of the distant rich cluster CI 0939+4713 [z — 0.41). We compare the optical and 
near-infrared morphologies of cluster members and find apparent small-scale optical structures within the 
galaxies which are absent in the near-infrared. We conclude that strong dust obscuration is a common 
feature in the late- type galaxies in distant clusters. We then concentrate on a sample of ten faint radio 
galaxies lying within our NICMOS field and selected from a very deep 1.4-GIIz VLA map of the cluster 
with a 1 (T flux limit of 9 /iJy. Using published data we focus on the spectral properties of the eight radio- 
selected cluster members and show that these comprise a large fraction of the post-starburst population 
in the cluster. The simplest interpretation of the radio emission from these galaxies is that they are 
currently forming massive stars, contradicting their classification as post-starburst systems based on the 
optical spectra. We suggest that this star formation is hidden from view in the optical by the same 
obscuring dust which is apparent in our comparison on the optical and near-infrared morphologies of 
these galaxies. We caution that even in the restframe optical the effects of dust cannot be ignored when 
comparing samples of distant galaxies to low-redshift systems, particularly if dust is as prevelant in 
distant galaxies as appears to be the case in our study. 

Subject headings: galaxies: evolution — galaxies: clusters: individual CI 0939+4713 (A 851) — 
galaxies: starburst — radio continuum: galaxies — infrared: galaxies 



I. INTRODUCTION 

The characteristics of the galaxy populations in the 
cores of rich clusters appear to vary strongly as a function 
of redshift, with an increasing fraction of blue galaxies in 
clusters ai z > 0.2 (the Butcher-Oemler effect, Butcher & 
Oemler 1984). The search for the processes responsible for 
these rapid changes has become a vigorous field of research 
and at least for the blue cluster galaxies a broad concen- 
sus has been reached about the nature of this population. 
The bulk are star-forming cluster members, they cover a 
wide range of luminosities, from a few L* downwards, and 
have emission line strengths which indicate moderate to 
high star formation rates (Abraham et al. 1996; Balogh 
et al. 1999; Barger et al. 1996; Couch & Sharpies 1987; 
Dressier & Gunn 1992; Dressier et al. 1999, D99; Fabri- 
cant et al. 1994; Fisher et al. 1998; Lavery & Henry 1988). 
The optical morphologies of these galaxies were investi- 
gated firstly with high-resolution ground-based (Lavery et 
al. 1992) and more recently with Hubble Space Telescope 
(HST) imaging (Dressier et al. 1994a; Couch et al. 1994, 
1998; Oemler et al. 1997; Small et al. 1997, S97; Fabri- 
cant et al. 1999). These studies have shown that most 
of the blue star-forming galaxies have strong disk com- 
ponents and a large fraction also appear to be disturbed, 
with structures suggesting mergers or tidal interactions. 
Thus the increasing proportion of blue galaxies within dis- 
tant clusters is associated with an increasing number of 



actively star-forming disk galaxies in these environments. 
The large fraction of apparently disturbed galaxies sug- 
gests that galaxy-galaxy interactions or the effects of the 
cluster tidal field on these disk galaxies could be respon- 
sible for the Butcher-Oemler effect (Moore, Lake & Katz 
1998). 

The spectroscopic studies of distant clusters uncovered 
another population of luminous galaxies which is not seen 
in similar numbers in local clusters. These are galaxies 
with strong Balmer absorption (typically measured us- 
ing the H(5 A4104A line and coming predominantly from A 
stars) and no [On] A3727A emission. The large population 
of A stars is a sign that the galaxy was actively forming 
massive stars in the recent past (< ICyr), while the lack 
of emission lines suggests that this star formation has now 
ceased. Together the spectral properties are interpreted 
as showing that the galaxy is in a post-starburst phase. 
Almost 20% of the luminous galaxies within the cores of 
distant clusters fall in this post-starburst class (e.g. D99, 
although see Balogh et al. 1999) and it has been suggested 
that this high fraction indicates that almost all cluster 
galaxies pass through this phase (Barger et al. 1996). 

There is an increasing realisation in the community of 
the importance of dust obscuration in defining the appar- 
ent properties of galaxies in the near-UV and optical, par- 
ticularly at high redshifts. This issue has been highlighted 
by the disagreement over the form of the star formation 
history of the Universe as estimated from the variation in 
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the ultraviolet luminosity density with redshift (e.g. Lilly 
et al. 1996; Steidcl et al. 1999) compared to star forma- 
tion indicators at longer wavelengths such as Ha emission 
(Yan et al. 1999) or reradiated star light detected in the 
sub- millimeter (e.g. Blain et al. 1999). The longer wave- 
length tracers tend to find higher star formation densities, 
a result which has been attributed to dust absorption in 
the UV and optical (Calzetti & Heckman 1999; Meurer et 
al. 1997; Pettini et al. 1998). 

The issue of the effects of dust on the evolutionary cy- 
cle associated with the optical Butcher-Oemler effect has 
been disregarded by most investigations to date. How- 
ever, the recent analysis of the spectral catalog for the 
MORPHS sample (D99) by Poggianti et al. (1999, P99) 
has uncovered evidence for dust obscuration in at least 
one spectral class - e(a) galaxies. These galaxies have 
detectable [On] emission and relatively strong Balmer 
absorption (EW(H(5) > 4). These spectral features, 
along with other properties of the galaxies and of similar 
spectrally-classified local systems (P99; Poggianti & Wu 
1999), were interpreted by P99 as probably arising from a 
dust-obscured star burst, with the [On] line strength sup- 
pressed by dust absorption. The e(a) class is naturally 
identified as the progenitor of the large population of post- 
starburst galaxies in the clusters. 

To obtain a complete view of the star formation prop- 
erties of cluster galaxies free from the effects of dust ob- 
scuration we need to complement the optical studies with 
observations at longer wavelengths, from the near-infrared 
to the sub-mm and radio. In particular, observations in the 
sub-mm and far-infrared can trace the amount of star light 
absorbed by dust and reradiated at longer wavelengths. 
For the most extreme starburst galaxies in the local Uni- 
verse this reradiated emission dominates the bolometric 
luminosity of the galaxies (e.g. Sanders & Mirabel 1996). 
Unfortunately, current sub-mm and far-infrared surveys 
lack the sensitivity to study all but the most extreme star- 
burst galaxies (M ^ lOOAfQ yr^-'^) at moderate and high 
redshift. However, at radio wavelengths, the synchotron 
radiation from electrons accelerated by supernovae leads 
to a tight correlation between the radio and far-infrared 
fluxes of star- forming galaxies across a wide range in lumi- 
nosity (Condon 1992). This means that in the absence of 
a radio-loud active nucleus the radio flux of a star-forming 
galaxy can be used to estimate its current massive star 
formation rate. 

By employing the huge collecting area of the VLA we 
can undertake a sensitive radio survey for star formation 
in obscured galaxies down to limits of M ~ 5-10Mq yr~^ 
at 2 ~ 0.5. To this end the VLA has been used in an 
ambitious survey to determine the evolution of the radio 
populations in a sample of very-rich clusters out to z = 0.4 
(Morrison et al. 1999). The survey contains 14 clusters at 
z = 0.10-0.18, 12 more between z = 0.20-0.25 and a fur- 
ther four at higher redshifts, z = 0.38-0.41 (the radio ob- 
servations discussed in this paper are of the most distant 
cluster in this survey). Morrison et al. (1999) find that 
the population of radio galaxies with luminosities above 
2 X 10^^ WHz"-'^ in the most extreme clusters has increased 
by a factor of ~ 5-6 out to z = 0.4. They conclude that the 
radio galaxy population in very rich clusters has evolved 
significantly over the last 5 Gyrs. 

In this paper we report on high resolution near-infrared 



imaging from HST and sensitive radio observations with 
the VLA of a small field in the core of the distant cluster 
C10939-h4713 (A 851, z = 0.41) for which high-resolution 
optical HST images and extensive archival spectroscopic 
observations are available. We combine these datasets to 
study the characteristics of galaxies in distant clusters and 
to search for the signatures of dust obscuration on their 
apparent properties. 

In §2 we describe the HST optical and near-infrared 
imaging as well as the deep VLA 1.4-GHz map of this 
field. In §3 we then present our results on the properties 
of the radio-selected galaxies within our HST field drawing 
on published spectroscopy. We discuss these in §4 and give 
our main conclusions in §5. Throughout this paper we use 
Ho = 50kms~"'^ Mpc^-'^ and a ilo = 1, A = cosmology. 
In this geometry 1" corresponds to 6.5 kpc at z = 0.41. 

2. OBSERVATIONS 
2.1. NICMOS Imuging 

The near-infrared NICMOS observations of 
CI 0939+4713 discussed here were obtained with Cam- 
era 3 during the special campaign in January 1998. Due 
to the distortion of the NICMOS dewar the wide-field 
Camera 3 (field-of-view 51.2" square, sampling 0.20" per 
pixel) was not par-focal with the other instruments on- 
board HST and to allow the diffraction-limited operation 
of this camera a special campaign was organised during 
which the main telescope optics were refocused on the 
Camera 3 focal plane. 

The observations comprise a mosaic of four pointings 
through the F160W filter covering a combined field of 
around 100" x 100". Each pointing consists of four indi- 
vidual exposures of 1215 s, each spatially offset by 1.1" on 
a square grid to remove bad pixels and reduce the effects 
of flatfield variations. The total integration time per pixel 
for the mosaic is thus 4.86 ks. The data were reduced in a 
standard manner using the NICMOS data pipeline (CAL- 
nica) and the most recent calibration files taken from the 
archive. Any remaining features in the background were 
removed with the unpedestal software kindly provided 
by Dr. R. van der Marel which worked well, and the final 
mosaicing was achieved using IMCOMBINE in IRAF. The 
mosaic was photometrically calibrated to give magnitudes 
on the i?i6o Vega system. This image provides restframe 
J-band morphologies at a resolution of around 1.3 kpc for 
cluster members (Fig. 1). 

To create a complete catalog of objects detected within 
the NICMOS mosaic we analysed the frame using the 
SExtractor image analysis package (Bertin & Arnoults 
1996). We adopted detection criteria of a minimum area 
of > 0.4 sq. arcsec (10 pixels) above an /xieo = 24.0 mag 
arcsec"^ isophotc and using these we detect 303 sources 
brighter than the approximate 3 a point-source flux limit 
of iJieo = 23.5. 

2.2. Radio Data 

The VLA observations of CI 0939-^4713 were obtained in 
B configuration at 1.4 GHz between 1996 January 06-08. 
The total integration time of the map is 60.7ks at an effec- 
tive resolution of 5.0" (33 kpc at z = 0.41). The map was 
cleaned and analysed using aips. Details of this complex 
procedure are given in Morrison et al. (1999). The map 
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reaches a 1 cr noise level of 9/iJy beam~^. We detect 10 
sources above 27 //Jy (3 a) within the field defined by the 
overlap region of our NIC MO S F160W mosaic and Wide 
Field and Planetary Camera 2 {WFPC2) F702W expo- 
sures (see Fig. 1). All ten sources have bright optical/near- 
infrared counterparts in our HST images. The properties 
of these ten sources are discussed in §3. The positions, 
fluxes and restframe radio power (using the redshifts from 
Table 2 and assuming an S oc spectrum) are re- 

ported in Table 1. Unfortimately due to the modest reso- 
lution of the map only two sources, #122 and #296, ap- 
pear to be resolved and then only marginally. The radio 
and optical/near-IR images are aligned to an rms accuracy 
of 0.4" using the 8 brightest radio sources in the NICMOS 
fleld. 

The flux limit of the radio sample corresponds to a radio 
power of just 2 x 10^^ W Hz~^ at the redshift of the cluster. 
Comparing this limit to those achieved by radio surveys of 
local clusters we note that Gavazzi & Boselli (1999) detect 
galaxies as faint as 3 x 10^" WHz^^ in their analysis of the 
radio luminosity function in the Virgo cluster (a factor of 
roughly lOOx closer to us than CI 0939-^4713) using the 
NVSS 1.4 GHz survey (Condon et al. 1998). 

2.3. WFPC2 Imaging 

These observations were obtained as part of the Early 

Release Observations during the Science Verification 
Phase after the installation of the WFPC2 on board HST. 
The observations comprise a total of 21.0 ks of integra- 
tion through the F702W filter (Dressier et al. 1994b). The 
exposures were dithered by integer pixels to allow the re- 
moval of defects and cosmic ray events. The position of the 
corrective optics within WFPC2, to correct for the aberra- 
tion of the primary mirror, was not optimal during these 
observations and combined with a higher than ideal oper- 
ating temperature for the CCDs these images are not as 
cosmetically clean as later exposures with WFPC2. Never- 
theless, they provide deep (R ^ 26) imaging of the core of 
the cluster at 0.1-0.2" resolution in the restframe B-band. 

The reduction and analysis of the F702W exposure of 
this field is described in more detail in Smail et al. (1997), 
who provide a complete catalog of objects detected in the 
fleld with photometry calibrated to the Rro2 passband de- 
scribed by Holtzmann et al. (1995). SmaU et al. (1997) 
also discuss visually classified morphologies on the revised 
Hubble scheme for the brighter galaxies within this field 
down to i?702 = 23.5. There total of 123 galaxies with 
morphologies from S97 within the joint NICMOS/ WFPC2 
fleld. All ten of the radio sources selected from the VLA 
map are included in the S97 catalog and we list the Hubble 
type for each galaxy in Table 2. We also give the Distur- 
bance Index (D), which is a measure of how disturbed the 
galaxy seemed compared to the typical appearence of a 
galaxy with its Hubble type locally, where £> = is "nor- 
mal" and D = 4 is highly disturbed (S97). In addition 
to these objective measures, a subjective interpretation of 
the source of the disturbance was also given by S97 and 
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is listed in the table (e.g. Chaotic, Tidal interaction, etc). 
We note in particular that a galaxy would be classed as 
"Merger" if two or more close nuclei were seen in a com- 
mon envelope. These visual estimates were shown to cor- 
relate well with the machine-based asymmetry measure- 
ments (S97) and have been used to attempt to isolate and 
study the role of dynamical interactions in triggering the 
Butcher-Oemler effect (e.g. D99). 
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Fig. 2. The (-R702 — ^^^leo) — ^leo color-magnitude distribution for the 
sources lying in the joint WFPC2/ NICMOS field. The (i?702 - Hieo) 
color is measured witliin an aperture of diameter 2.5 times the FWHM 
of the galaxy in the F160W image. The large filled points denote those 
galaxies detected in the VLA map, the circled points are other galaxies 
with spectroscopic identifications in D99 and the small points are the 
remaining sources selected from the F160W catalog. Note the sequence 
of red galaxies associated with early-type cluster members which extends 
down to HiQo ~ 23. 

To investigate the optical-infrared colors of the galaxies 

in our sample we have aligned and resampled the WFPC2 
image to the reference frame of the NICMOS exposure. A 
comparison of stars and compact galaxies in the resam- 
pled F702W exposure and the F160W image shows that 
they have very similar proflles in the two bands indicat- 
ing that the resolution is dictated by the pixel sampling 
of the images. This suggests that PSF variations should 
not strongly influence our analysis. To confirm this we 
have convolved the WFPC2 and NICMOS images with 
a model PSF for the other instrument generated using 
TINYTIM (Krist & Hook 1999). The colors determined 
from these convolved images show no systematic differ- 
ences from those measured from the unsmoothed images 
and hence we use the latter in the following analysis. 

To measure colors from our F702W and F160W expo- 
sures we have used the IRAF task phot and apertures 
matched to the near-infrared extent of the galaxies. The 
diameter of the photometry apertures were taken as 2.5 
times the FWHM of each source from the F160W cata- 
log (dphot in Table 2). The color-magnitude diagram for 
the fleld is shown in Fig. 2, with the radio sources and 
spectroscopically identifled galaxies marked. We give to- 
tal HiQQ magnitudes and aperture (i?702 — -ffieo) colors for 
the ten radio sources in Table 27 We show the F160W 
and F702W images of this sample in Fig. 3. 

To study the color variations within the galaxies we have 



'^For interest we list here the (R702 — ^^l6o) colors of the various z > 3 lensed galaxies identified by Trager et al. (1997). These galaxies are 
all relatively blue compared to the general field (Fig. 2), (^702 - i^ieo) = 1-25 ± 0.10 for DG#334 aX z = 3.34, (i?702 - -ffieo) = 0.78 it 0.15 
for the z = 3.97 source seen as P1/P2/P3 and {Rro2 — Hiqo) = 1.43 ± 0.20 for the z = 3.98 arclet AO, although this may suffer contamination 
from the nearby bright cluster galaxy (#426, Fig. 3). A search of the NICMOS and WFPC2 images to identify galaxies with extreme colors 
turned up a very red disk galaxy near the center of the field: #333 (Fig. 1), this is the reddest, bright object in the field with Hiso = 19.53 
and (R702 - Hieo) = 5.4 ± 0.2 (c.f. Treu et al. 1999). 
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Table 1 

Properties of the Radio Galaxy Sample 



Source 


R.A. 


Dec. 


5'l.4 


Pi. 4 


Comments 


ID 


(,12000) 




(10^^ WHz-1) 




36 


09 42 58.16 


+46 59 53.4 


79 ± 9 


5.9 




122 


09 42 57.62 


+46 59 45.3 


204 ± 9 


15.2 


Resolved? 


224 


09 42 55.80 


+46 59 40.4 


55 ± 9 


4.0 




230 


09 43 02.32 


+46 59 28.7 


52 ± 9 


3.8 


Pi. 4 assumes z = 0.41 


296 


09 43 02.81 


+46 59 24.2 


644 ± 9 


47.6 


Resolved? 


399 


09 42 57.84 


+46 59 13.3 


43 ± 9 


3.1 




426 


09 42 56.21 


+46 59 12.0 


32 ±9 


2.4 




536 


09 43 00.03 


+46 58 53.4 


196 ±9 


14.5 




610 


09 42 57.34 


+46 58 50.4 


158 ±9 


11.7 




640 


09 42 55.08 


+46 58 41.7 


57 ±9 


2.6 





constructed F702W— F160W color images for each galaxy 
in our spectroscopic sample. We show the color images for 
the relevant radio-detected galaxies in Fig. 3 to illustrate 
the internal color variations within these galaxies. A com- 
parison sample of confirmed cluster members which are 
not detected in the radio map is shown in Fig. 4. We note 
that PSF-convolved versions of these images are qualita- 
tively similar, in particular the blue cores at the centers 
of #296, #399 and #426 do not appear to result from 
differences in the PSF between the two passbands. 

2.4. Spectral Information 

Spectroscopic observations of galaxies within our field 
are available from the MORPHS spectral catalog (D99). 
The targets for these observations were selected from the 
F702W exposure and D99 amassed spectra for 28 galaxies 
within the NICMOS field, comprising 22 cluster members 
and 6 field galaxies. These spectra typically cover the 
spectral region from 3000-6000A in the restframe, includ- 
ing any [Oii]A3727 emission but in most cases not Ha. 
The galaxies have all been classified using the spectral 
classification scheme of D99 as follows: passive galaxies 
(k); post-starburst (k+a/a+k), these show no [On] and 
strong Balmer absorption indicative of a cessation of star 
formation within the last Gyr or so; actively star forming 
(c(c)/c(b)/c(a)), e(c) have spectral properties similar to 
local spiral galaxies, c(b) have stronger [On] and arc inter- 
preted as starburst galaxies, e(a) have moderate [On] but 
stronger Balmer lines than e(c) galaxies. The e(a) class 
has been interpreted as dusty starbursts (P99) where the 
[On] is suppressed by dust within the galaxy (supported 
by the observations of anomalously high Hq:/[Oii] ratios 
for these galaxies, see also §1). The distribution of clus- 
ter/field galaxies in the different spectroscopic classes is: 
k (7/2); k-ha (6/0); a+k (2/0); c(a) (5/1); e(b) (1/0); e(c) 
(2/1). The spectrum of one field galaxy has too low signal- 
to-noisc to allow a reliable classification. 

For the ten radio sources in the field, spectra and spec- 
tral classifications are available from D99 for nine galaxies 
(galaxy #230 lacks a spectroscopic ID) and wc list these in 
Table 2. Of these nine galaxies, eight are cluster members 
and one is a field galaxy (the Sab galaxy #640). 

3. RESULTS 



With a small sample of galaxies our conclusions are by 
necessity fairly qualitative in nature and this paper should 
be seen as the prelude to the analysis of a larger sample of 
radio-selected galaxies in distant clusters by Morrison et 
al. (1999). Nevertheless, we will use our observations, and 
particularly the high resolution near-infrared imaging with 
NICMOS, to illustrate the general characteristics of radio 
galaxies in distant clusters and to emphasise the effects of 
dust in galaxies on their perceived properties. 

3.1. Radio-selected Sample 

We have detected ten galaxies above the 3(t flux limit 
{Si 4 > 27/iJy) of the VLA map which lie within the 
joint NICMOS/ WFPC2 field. The spectroscopic infor- 
mation available from D99 allows us to identify eight of 
these galaxies as cluster members (or 35% of the 22 clus- 
ter galaxies known within this field) and a further source 
as a foreground field galaxy (12% of the spectroscopically- 
confirmed field population within our image). Our field 
size corresponds to a region ~ 650 x 650 kpc in extent at 
the cluster redshift and the equivalent radio power limit is 
2 X 10^^ WHz^^. A comparably deep survey of a similar 
region in the core of a rich local cluster would typically 
detect ^ 1 galaxy (e.g. Fig. 5), compared to the eight de- 
tected here. This suggests a substantial increase in the 
rate of occurence of radio sources in high density environ- 
ments at moderate redshifts (Morrison et al. 1999). 

The galaxies included in the radio-selected sample of 
cluster members are typically the brighter galaxies in our 
field (Fig. 2), including the brightest four galaxies (in the 
restframe i?-band). The whole sample spans a range of 
magnitudes i?702 = 18.3-20.8 {Lb ^ A-QAL%). The col- 
ors of the fainter radio-detected galaxies are slightly bluer 
than equivalently luminous, undetected cluster members, 
with the brighter radio-emitting galaxies having compara- 
ble (red) colors to the undetected population (Fig. 2). The 
fact that we only detect the optically brightest galaxies in 
the radio map may suggest that the majority of cluster 
galaxies would be detected if we could increase the sensi- 
tivity of the radio observations by a factor of a few (see 
Fig. 5).8 

In our analysis we will assume that the local relation- 
ship between star formation rate and radio emission (e.g. 
Condon 1992) holds for spiral galaxies detected in distant 
clusters. The validity of the local relationship at least 

*We note that the radio population has similar dynamics to those cluster members undetected in the radio map: we calculate a restframe 
velocity dispcrion of (7 = (960 ± 160) kms~^ for the eight radio-detected members compared to ct = (1110 ± 140) kms~^ for the 12 undetected 
members lying within the same field. 
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Table 2 

Optical/Near-infrared Properties of the Radio Sample 



Source 


z 


.^160 


(.R702 — ^^i6n) 


'^phot 


Morph. 


Spectral 


EW(H<5) 


D 


Comments 


ID 








(") 


Type 


Type 


(A) 






36 


0.4119 


18.20 


1.46 ± 0.10 


5.6 


Sd 


e(a) 


10.2 


2 


EW([Oii])= -6.9A; Merger 


122 


0.4125 


17.96 


2.30 ±0.10 


2.7 


Sd 


a+k 


11.6 


3 


Merger, H(5 measurement uncertain. 


224 


0.4076 


15.85 


2.33 ± 0.04 


5.8 


Scd? 


k+a 


6.1 


2 


Merger 


230 




18.80 


1.77 ± 0.14 


2.2 


Irr 






2 


no spectrum in D99 


296 


0.4014 


16.87 


2.49 ± 0.07 


2.5 


E 


k+a 


4.6 







399 


0.4109 


16.14 


2.50 ± 0.05 


2.5 


E 


k 


0.0 







426 


0.4037 


16.18 


2.55 ± 0.05 


3.9 


E 


k 


0.0 







536 


0.4111 


17.49 


2.29 ± 0.08 


4.0 


So 


a+k 


8.6 


3 


Tidal? 


610 


0.4007 


15.84 


2.50 ±0.05 


3.4 


Sa/SO 


k+a 


4.2 


1 


Merger? 


640 


0.3324 


16.27 


3.21 ± 0.04 


8.3 


Sab 


k: 


0.0 


1 


Field galaxy, low s/n spectrum — possible e(a)? 



for the distant field appears to be supported by recent 
observations (Georgakakis et al. 1999). However, we cau- 
tion tliat tlic exact star formation rates derived from this 
analysis may have substantial systematic uncertainties. In 
particular the claim that the radio to far-infrared correla- 
tion varies locally between galaxies in rich clusters and 
those in the field (Andersen & Owen 1995) may indicate 
an environmental influence in the conversion from radio lu- 
minosity to star formation rate (SFR) resulting from the 
compression of the galactic magnetic field by the galaxy's 
motion through the intraclustcr medium (ICM). Without 
a fuller understanding of the state of the ICM in distant 
clusters it is difficult to quantify the extent of this bias 
on the radio-derived SFRs of distant cluster galaxies. We 
have chosen therefore to quote the SFRs determined us- 
ing the local field relationship and the correction which 
should be applied if the radio population in distant clus- 
ters shows the same behaviour as is claimed for galaxies 
in rich clusters locally. 

3.2. Optical Morphologies 

Turning to the restframe i3-band morphologies of the ra- 
dio sources (Table 2 & Fig. 3) we see that they fall roughly 
equally into two broad classes: optically- luminous early- 
type galaxies (E and SO/Sa) and mid- to late- type spiral 
galaxies (Sc-Sd/Irr). 

In their analysis of galaxies in the Virgo cluster detected 
in the NVSS radio survey, Gavazzi & Boselli (1999) found 
that the most powerful radio sources in spirals were also 
in mid-type, Sb Sc, systems. To compare the optical and 
radio properties of the galaxies we plot the 1.4-GHz rest- 
frame power versus the optical B-band luminosity for all 
spectroscopically-confirmed members in Fig. 5, separated 
into early- and late- type. The early- type galaxies show 
a wider range in their Pi ^/Mb values than the spirals, 
but span a similar region of the Pi. 4-Mb plane as radio- 
selected E/SO galaxies in the Virgo cluster (Fig. 5). This 
is as expected if the radio emission from the early-type 
galaxies is coming from a central AGN. In contrast the 
radio emission from the spirals is (by analogy to local sys- 
tems) most likely synchotron emission associated with ac- 
tive star formation and leading to a correlation between 
the radio power and the optical luminosity of the galaxy. 
We note that the population of radio-emitting spiral galax- 
ies identified by Gavazzi & Boselli (1999) in Virgo have 
radio powers 5x lower at a fixed B-band luminosity than 
the galaxies we are detecting in CI 0939-1-4713 (Fig. 5). 
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Fig. 5. The Pi,4-Ms distribution for those galaxies detected in the 
VLA map of CI 0939+4713 compared to the equivalent distribution for 
the Virgo cluster (Gavazzi &: Boselli 1999). The top panel compares the 
distributions for the morphologically-classified spiral and irregular galax- 
ies in both clusters, while the lower panel shows the equivalent compar- 
ison lor elliptical and SO members. We adopt an apparent magnitude of 
i?702 — 19.7 for an Mb — —21.1 (L^) galaxy in the i?702-band observed 
at 2: — 0.41 (assuming an SED typical of a present-day Sbc galaxy). The 
dotted lines show the radio flux limit of our sample in CI 0939+4713. 
We indicate the spectroscopically-confirmed member galaxies which are 
undetected in the VLA map as upper limits. 

The association of around 60% of the bright radio- 
selected sample with spiral galaxies suggests that the evo- 
lution in the radio populations within clusters runs in par- 
allel to the increasing fraction of blue star- forming galaxies 
found in these regions at 2; ^ 0.2, the Butcher-Oemler ef- 
fect. The rise in the fraction of bright spiral galaxies in 
clusters such as CI 0939-1-4713 is a factor of ^ 5 out to 
z = 0.4 (Dressier ct al. 1997). Even without the increased 
radio luminosity of the brightest cluster spiral galaxies, 
this increase in the proportion of spirals is enough to ex- 
plain the evolution in the number of radio sources seen 
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in the distant clusters (a factor of ^ 5 6, Morrison ct al. 
1999). Tlic increase in the spiral fraction in the clusters is 
believed to result from infall from the surrounding field as- 
sociated with the growth of the clusters. A rise in the radio 
activity of moderate redshift field galaxies has also been 
observed and attributed to star formation in these galax- 
ies (Georgakakis et al. 1999). However, to fully understand 
the evolution of the radio population in the cluster we also 
need to identify the cause of the increased radio luminosi- 
ties of the bright spiral galaxies. We note that our use of 
a local cluster sample in the comparison should minimise 
the possible environmental effects which could explain this 
apparent brightening, and moreover the increase in radio 
luminosities is substantially above that claimed to result 
from environmental effects within local clusters (§3.1). 

3.3. Near-infrared Morphologies 

We now discuss the internal structure of the cluster 
galaxies in our HST images (Fig. 3 and 4) and the im- 
portant insights which our NICMOS near-infrared imag- 
ing has provided into the nature of these systems. Start- 
ing with the early-type galaxies we see relatively uni- 
form optical-near-infrared colors (as shown by their lack 
of structure in the F702W— F160W images in Fig. 3 and 
4). A few contain blue cores, which may be associated 
with central AGN, although the details of these features 
are sensitive to the exact form of the PSF in the F702W 
and F160W images and so we do not discuss them further 
here. 

A comparison of the F702W and F160W morphologies 

for the late-type radio sources in Fig. 3 shows that in these 
systems there are substantial differences in the appearence 
of these galaxies between the restframc; B- and J-bands. 
The near-infrared morphologies are considerably more re- 
lax;ed and symmetrical. However, the lack of a morpholog- 
ical classification scheme for local galaxies based on near- 
infrared imaging means that is not useful to attempt to 
reclassify the galaxies on the basis of the F160W images. 

More interestingly, the F702W-F160W images in Fig. 3 
show that in a large number of cases the irregular struc- 
tures seen in the rest-frame B-band within the central re- 
gions of the late-type galaxies are replaced by a much more 
regular appearence in the near-infrared and so must be due 
to very red structures (visible as white regions in Fig. 3). 
These features are substantially redder than the outskirts 
of the galaxies (up to A(i?— _ff) ^ 1.3 1.5) and could repre- 
sent differences in the stellar populations within the galax- 
ies or simply dust obscuration. To distinguish between 
these alternatives we concentrate on the large fraction of 
galaxies whose integrated spectra show post-starburst stel- 
lar populations. The lack of emission linc;s in the; optical 
spectra indicates that there is no visible star-formation in 
the outskirts of these galaxies, which also show the red 
optical-infrared colors typical of an evolved stellar popu- 
lation (see the integrated colors in Fig. 2). In this sit- 
uation it is very difficult to make the stellar population 
in the central regions of the galaxies substantially redder 
than the outskirts without invoking a contrived scenario. 
In contrast, dust obscuration can easily produce both the 
apparently red colors and its structured distribution. We 
propose therefore that large quantities of dust are proba- 
bly a common feature of these galaxies. 

The small scale features seen in the optical images of 



these galaxies which wc claim are due to dust obscuration 
have been previously interpreted as arising from dynami- 
cal disturbance of these galaxies, suggesting a dynamical 
origin for their star formation activity (S97). To under- 
stand the consequences of this change in interpretation 
of these galaxies we focus on those galaxies classified in 
the optical as showing signs of morphological disturbance 
by S97. There are a total of eight spectroscopically con- 
firmed members in our sample with disturbance indices of 
D > 2 indicating abnormally disturbed morphologies, all 
of these are disk galaxies and they are typically classified as 
mergers or tidally disturbed, half of them are detected in 
our radio map. These galaxies exhibit substantially more 
regular morphologies in the near- infrared (Figs. 3 & 4), 
in particular the apparently double structures seen in the 
F702W morphologies of the radio sources #36, #122 and 
#536 all disappear in the F160W passband. 

We also see signs of dust obscuration in the disturbed 
galaxies which are not radio sources (#171, #369, #497 
and #622) , although arguably not as extreme as the radio- 
selected sample. In a few cases there are other signatures 
of interactions in the F702W images apart from double 
nuclei, etc, such as low-surface brightness tidal features 
or companion galaxies (e.g. #224 or #497), which indi- 
cate that these are true merging or interacting systems. 
But, in the majority of the galaxies their high disturbance 
index was based upon the apparent presence of multiple 
nuclei in the central regions of the galaxies (see §1). We 
suggest that in fact dust obscuration is the correct ex- 
planation for their irregular appearence in the restframe 
B-band which results in their classification as disturbed or 
merging. We estimate that of the eight confirmed mem- 
bers with D > 2 based on the F702W imaging, at least half 
(#36, #122, #171, #622) would have £) < 2 if classified 
from the F160W images, suggesting that the occurence of 
highly disturbed galaxies in the core of CI 0939-1-4713 may 
have been overestimated by up to a factor of two. 

We conclude from our comparison of the near-infrared 
and optical morphologies of a small sample of confirmed 
cluster members that the frequency of galaxies showing 
small-scale signatures of disturbance or interaction is prob- 
ably substantially less than has been suggested from opti- 
cal studies (e.g. S97) due to obscuration by dust. This 
effect acts in addition to the standard luminosity bias 
in optically-selected samples towards actively star- forming 
galaxies (the usual justification given for undertaking near- 
infrared surveys of distant clusters, e.g. Barger et al. 1996). 
The expectation of previous optically-based studies was 
that dust obscuration would not be a significant bias when 
comparing the distant morphological samples with more 
local examples as they were comparing effective restframe 
i?-band observations in both cases (S97). However, this 
expectation seems to have been misleading due to the 
wider prevelance of dust in cluster galaxies at z ~ 0.4, 
perhaps associated with the increased activity in galaxies 
at these epochs. 

3.4. Spectral Properties 

Looking now at the spectral properties of the radio- 
selected cluster members we see something striking: five 
of the galaxies have spectral features which class them 
as post-starburst (a-|-k/k-|-a) systems using the scheme of 
P99 (Table 2). Moreover, these five radio sources show 
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the strongest Balmer absorption lines of the eight post- 
starburst galaxies in our cluster sample (they are also the 
brightest five post-starburst galaxies). There are no post- 
starburst field galaxies in our sample. 

The detection of post-starburst galaxies in the radio 
waveband is surprising. These galaxies are selected to 
have no detectable [On] emission and hence are expected 
to have no on-going massive star formation (observations 
of the spectral region including any Ha emission in these 
galaxies would make this statement more robust). In con- 
trast, the radio emission is most easily explained as arising 
from massive star formation (although some contribution 
from an obscured radio-loud AGN cannot be ruled out 
with our present data). The expected lifetime of the radio 
emission after the starburst is ^0.1 Gyr, thus if they are 
post-starburst systems these galaxies must all be seen just 
after the cessation of their star formation. This seems very 
unlikely as the length of time in which the post-starburst 
signature (strong B.5 absorption) is visible is ~ 1 Gyr and 
we would thus expect to see an order of magnitude more 
radio-quiet k-|-a galaxies, where the radio emission had 
decayed but the Balmer absorption is still visible. This 
population is not observed. Therefore, we conclude that 
massive stars are currently being formed in these galaxies. 

Where is the site of the active star formation in the 
a-hk/k-ha cluster galaxies? The optical spectra typically 
sample the whole of the bulge and disk of the galaxy (see 
D99) and thus it is unlikely that we are missing emission 
from the outskirts of the galaxy. We suggest instead that 
the star formation is occuring in the central regions of 
these galaxies but is hidden from view by the dust which 
we see there in our F702W— F160W images. This pos- 
sibility was raised in P99 who stated that the most ex- 
treme e(a) galaxies might appear as post-starburst sys- 
tems {a.+k/k+a.) due to obscuration of the emission re- 
gions. The reddened regions lie in the central 5-10 kpc of 
the galaxies (Fig. 3 & 4) and have optical-near-infrared 
colors which arc typically A(i?702 — -ffioo) ~ 1 magnitude 
redder than the outskirts of the galaxies. This would in- 
dicate optical extinctions of at least ~ 2 mag to the 
population of disk stars (let alone the active sites of star 
formation) although this estimate is highly uncertain due 
to resolution effects. 

Four of the five radio-selected post-starburst galaxies 
are disk galaxies and we can estimate their star forma- 
tion rates (SFR) from their radio power using the calibra- 
tion of Condon (1992). We find that these galaxies are 
typically forming massive stars {M > 5Mq) at rates of 
~ 10 M© yr-\ with the highest SFRs being 30 Mq yr"^ 
for #122 and #536. This assumes no boosting of the ra- 
dio emission due to compression of the galactic magnetic 
fields by the cluster ICM, which could amount to a fac- 
tor of ~ 2 (Andersen & Owen 1995). Combining these 
SFRs with the B-band luminosities of the galaxies we can 
predict their [On] fluxes if the star formation regions are 
unobscured using the relations given in Kennicutt (1992). 
We can then roughly estimate the amount of extinction 
needed for this line to be undetected in the optical spec- 
tra, EW([Oil])< 5A. We find that the star forming regions 
would need to suffer As ^ 3 mag of extinction for the [On] 
line to be undetected in the D99 spectra {Ab ^ 2 mag if 
we assume the radio fluxes are artificially raised). This 
amount of obscuration is not unreasonable given the ex- 



tent of the dust reddening seen in the F702— F160W im- 
ages (Figs. 3 & 4). We note that if we assume; this dust 
is cold, Td ~ 40 K, then the predicted sub-mm fiuxes of 
these galaxies will be only Ss5o ~ 0.5 mJy at 850 /um. This 
is well below the confusion limit of current sub-mm instru- 
mentation in blank fields. 

There axe also examples of proposed dusty galaxies 
within our sample which are not detected in the radio map. 
In particular three of the cluster galaxies have e(a) spec- 
tral features and evidence for dust in their F702— F160W 
images (Fig. 4, #177, #369 and #622). These galaxies are 
optically fainter than the e(a) galaxy which is detected in 
the radio, #36, and it is possible that their lower overall 
luminosities mean that their radio emission is below the 
limit of the VLA map. Interestingly, although the e(a) 
spectral class was interpreted as dusty starbursts in P99, 
the radio map places upper limits on the current rate of 
massive star formation in the undetected e(a) galaxies of 
M ^ 5Mq yr~^, with similar limits applying to the galax- 
ies with e(c) spectra in the field (e.g. #165 and #497 in 
Fig. 4). Taking the possible boosting of the radio emission 
due to compression of the galactic magnetic fields by the 
cluster ICM (§3.1) into account these SFR limits could be 
dropped by a factor of two or more. 

4. DISCUSSION 

We see a combination of post-starburst spectral fea- 
tures, radio emission and apparently large quantities of 
dust within the spiral galaxies in our sample. These char- 
acteristics can be most easily accomodated if they exist 
in spatially distinct regions of the galaxies, as suggested 
by our color maps. The radio emission would arise from 
massive; star formation in the dust obscured center of the 
galaxy, while the A stars trace the remanents of an earlier 
(^ 1 Gyr) strong star formation episode within the galactic 
disk. There is weak evidence in our sample for a correla- 
tion between US line strength and radio luminosity which 
would suggest that there is a causal link between the num- 
bers of A stars and the current (obscured) star formation 
rate, probably as a result of a single triggering mechanism. 

This triggering mechanism could be a cluster-related 
phenomenum, a view which is supported by the lack of 
obvious post-starburst systems in samples of faint radio 
sources in the field (Hammer et al. 1995; L.L. Cowie, priv. 
comm.). Thus, we may be seeing a more dynamic situ- 
ation including cluster-driven processes, whereby interac- 
tions between the intracluster medium and the cold molec- 
ular gas within the galactic disks accelerates the collapse 
and evaporation of the star forming regions. The effects 
of the ICM on the galactic disk are expected to be a func- 
tion of the loc;al density in the disk with lower density 
regions in the outskirts of the disks being most severely 
effected, while star formation continues almost uneffected 
in the inner regions of the galaxy. Thus environmental 
processes could result in a combination of post-starburst 
and starburst within individual cluster galaxies. Studies 
of the spatial distribution and kinematics of larger sam- 
ples of radio-selected galaxies may be able to discriminate 
between the different physical processes which could be 
operating (Morrison et al. 1999). 

Alternatively, perhaps this activity may be triggerc;d by 
a more general physical process, unrelated to the clus- 
ter environment. Indeed, similar behaviour is seen in 
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the models of interacting field galaxies published by Mi- 
hos & Hernquist (1996) (although we have argued that a 
substantial fraction of apparently disturbed cluster spiral 
galaxies probably arise from dust obscuration, rather than 
true interactions). For a range of models for tidal inter- 
actions between galaxies Mihos & Hernquist predict that 
spatially-extended enhanced star formation will first occur 
in the disk of the galaxy. This is followed on a timescale 
of a ~ 0.1-0.3 Gyr by funnelling of any remaining gas into 
the center of the galaxy where it reaches very high densi- 
ties and triggers an (obscured?) nuclear starburst. Again, 
more detailed investigation of the kinematics and internal 
dynamics of these galaxies may provide useful diagnostics 
to distinguish between the competing models. 

The picture which is appearing is one in which galax- 
ies which are visibly forming stars appear to have mod- 
est star formation rates, ^ 5Mq yr~-^ (M > 5Mq), while 
more vigorous star formation, ^ 10-30Mq yr"-"^ is on-going 
in highly obscured regions in the centers of what appear 
to be post-starburst galaxies. If the radio fluxes of these 
distant galaxies are increased by the same environmental 
effects which are claimed to boost the radio emission of 
local cluster galaxies and to the same extent (Andersen & 
Owen 1995) then the quoted SFRs should be reduced by 
a factor of ^ 2. Sensitive observations to search for Ha 
emission (or emission lines further into the near-infrared) 
will provide a stronger test of the star formation rates in 
these galaxies, the extent of obscmation within them and 
the degree to which their radio emission deviates from that 
seen in the distant field population due to environmental 
processes. Finally we note that by having post-starburst 
and starburst phases occuring in parallel within a single 
galaxy it is possible to reduce the high fraction of galax- 
ies (~ 100%, Barger et al. 1996) which must pass through 
this phase to explain the large numbers of post-starburst 
systems seen within distant clusters (Barger et al. 1996; 
P99). 

5. CONCLUSIONS 

In summary, we have highlighted the role that dust plays 

in obscuring our view of the morphologies and star forma- 
tion properties of distant galaxies. We have shown that se- 
lection at faint radio fluxes provides a powerful technique 
for identifying spectrally-classified post-starburst and star- 
burst galaxies crucial to our understanding of the evolu- 
tionary cycle of galaxies in clusters (e.g. P99). Using this 
new tool we can undertake statistical studies to search for 
the mechanism responsible for triggering and/or quench- 
ing of this activity (P99; Morrison et al. 1999). The main 
conclusions of our study are: 

• We present optical and near-infrared imaging with 
HST of a region in the core of the distant rich cluster 
C10939-F4713 at = 0.41. 

• We compare the high resolution near-infrared and opti- 
cal morphologies of a sample of 22 cluster members. This 
shows that while the broad morphological classifications 
are similar, the frequency of galaxies showing small-scale 
signatures of disturbance or interaction is substantially less 
when using the near-infrared imaging. We suggest that 
obscuration by dust has led optical studies to conclude 
that an anomalously high fraction of galaxies have suf- 
fered small scale disturbances. This may point to a wider 
prevelance of dust in cluster galaxies at ^; ~ 0.4, perhaps 



as a result of the increased star formation activity at these 
epochs. 

• Using a very deep VLA 1.4-GHz map of this area we se- 
lect a sample of ten radio galaxies in our field above a flux 

limit of 5*1.4 = 27/iJy, equivalent to a restframe 1.4-GHz 
power of 2 X 10^^ W Hz at the cluster redshift. The faint ra- 
dio sources are associated with two populations of cluster 
galaxies: luminous early-type galaxies, the radio emission 
from these galaxies originates from AGN; and dusty galax- 
ies with late- type spiral morphologies, the radio emission 
here is expected to trace massive star formation (although 
the exact relationship between the radio luminosity and 
the SFR is uncertain due to environmental influences). 

• The optical spectral characteristics of the late- type radio 
members are highly unusual - they appear to be post- 
starburst systems. We suggest that these galaxies are 
not post-starburst but in fact host highly obscured star- 
formation and star bursts. This is supported by our high 
resolution optical-near-infrared imaging which shows that 
the central regions of these galaxies are heavily dust en- 
shrouded. If they continued for several lOOMyrs, these 
dust enshrouded nuclear-starbursts could build a substan- 
tial stellar bulge in the remanent of these luminous disk 
galaxies, possibly contributing to the formation of the SO 
population in the clusters (Mihos & Hernquist 1996). 

• We emphasise that a lack of consideration of the effects 
of dust on the morphological and spectral properties of 
galaxies, even at the relatively modest redshifts discussed 
here, may lead to an incomplete understanding of the cy- 
cle of star formation in cluster galaxies. Further work on 
this subject is urgently required. 

• This paper is the prelude to a larger study of the HST 
morphologies and spectral properties of radio-selected 
galaxies in CI 0939-1-4713 and other high redshift clusters 
(Morrison et al. 1999). We expect that extensive studies 
with higher spatial resolution and wider wavelength cov- 
erage in the near future using the upgraded VLA and the 
next generation of millimeter interferometers will provide 
new insights into the prevelance and distribution of dust 
in distant galaxies. 
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Fig. 1. The VLA 1.4-GHz map overlayed as a contour plot on the HST NICMOS F160W mosaic. The radio sources are identified in large font using 
the numbering scheme of Small ct al. (1997). The galaxies labelled with the smaller font represent other spectroscopically identified galaxies (with the 
exception of :^333) from the catalog of Dressier et al. (1999), these are similarly numbered using the S97 identifications. The dashed line shows the 
overlap on this field of the WFPC2 F702W exposure from S97. The radio contours are at 12, 25, 50, 100 and 200 pjy per beam. North is top and east 
is left in this figure. 



Fig. 3. The groups of NICMOS F160W, F702W-F160W and WFPC2 F702W images for each of the ten galaxies selected from the VLA 1.4-GHz map 
of CI 0939+4713. The F160W and F702W images provide rcstframc J- and S-band views of the galaxies and these are shown with a linear intensity 
scale, while the color image is constructed from the difference of the log-scaled F160W and F702W images. In the F702W— F160W color image lighter 
regions represent redder areas of the galaxy, the typical color range is roughly -R702 — Hi&q — 2.4-3.2 (black to white). The panels are labelled with the 
catalog number, morphological type, spectral type and disturbance index (from S97 and D99). Each panel is 10" X 10" (equivalent to 65 kpc square at 
the cluster redshift) with north top and east left and the panels are ordered on morphology. Note that #640 lies in the extreme corner of the WFPC2 
field and #230 falls on a chip-boundary in the F702W exposure. All of the galaxies are confirmed cluster members with the exception of #640 which 
is a foreground field galaxy and #230 which has no spectroscopic identification. 



Fig. 4. A comparison sample of spectroscopically-confirmed cluster galaxies which show no detectable radio emission in the VLA map above a 1.4-GHz 
flux of Si. 4 = 27)LiJy. This figure includes all the spectrally-active members of this sample and again is ordered on morphology. The panels are labelled 
with the galaxy ID, morphology and disturbance index, D, from S97 and the spectral classification from D99. 
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